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1.  Overview 

The  research  supported  by  this  grant  was  focused  on  the  fundamental  physics  and 
applications  of  spin  polarized  atoms.  The  principal  investigator,  William  Happer,  is  very 
grateful  for  the  support  he  has  received  from  AFOSR  over  most  of  his  career  in  physics. 
At  the  end  of  the  academic  year,  2013-1014,  he  will  retire  from  teaching  at  Princeton 
University.  He  plans  to  continue  activities  in  physics  in  an  emeritus  status,  but  without 
the  graduate  students  and  post  docs  that  have  meant  so  much  to  him.  He  will  devote 
as  much  time  as  possible  to  the  physics  of  the  Earth’s  greenhouse  effect.  How  atoms 
and  molecules  emit  and  absorb  radiation  in  gases  with  pressure  ranges  and  temperature 
ranges  comparable  to  those  of  the  earth’s  atmosphere  has  been  a  central  focus  of  Happens 
AFOSR-supported  research,  and  will  be  an  great  advantage  for  work  on  greenhouse  physics. 
Of  particular  interest  is  how  collisions  of  CO2  molecules  with  the  nitrogen  and  oxygen 
molecules  of  the  air  affect  the  way  CO2  absorbs  and  emits  radiation.  Incorrect  modeling 
of  collisional  line  broadening  is  one  of  several  possible  reasons  that  climate  models  have  so 
greatly  exaggerated  the  warming  from  more  CO2.  There  has  been  basically  no  warming  at 
all  for  about  15  years,  contrary  to  model  predictions  of  about  0.3  C  temperature  increase 
for  the  increase  of  atmospheric  CO2  from  370  ppm  to  400  ppm.  Ultimately,  the  truth 
about  collisional  effects  will  have  to  be  resolved  by  incisive  experiments.  Theory  alone  is 
not  yet  sufficiently  quantitative. 

To  give  the  flavor  of  the  work  supported  by  this  grant,  in  the  next  section  we  will 
summarize  research  directions  of  particular  interest  to  the  principal  investigator,  with  a 
special  focus  on  the  more  theoretical  activities.  More  details  about  the  extensive  exper¬ 
imental  work  supported  by  the  grant  can  be  found  in  the  publications  listed  at  the  end 
of  this  report.  Also  listed  there  are  students  and  other  researchers  who  received  full  or 
partial  support  from  the  grant. 

2.  Some  Major  Accomplishments 

The  book,  Optically  Pumped  Atoms.  The  principal  investigator,  William  Hap¬ 
per,  and  two  former  students,  Thad  Walker,  now  a  professor  of  physics  at  the  University 
of  Wisconsin,  and  Dr.  Yuan- Yu  Jau  now  a  member  of  the  professional  staff  at  Sandia  Na¬ 
tional  Laboratory,  wrote  a  book  “Optically  Pumped  Atoms,”  published  by  Wiley  in  2010 
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[8]  which  summarizes  many  of  the  powerful  new  theoretical  methods  that  have  been  devel¬ 
oped  for  modeling  the  evolution  of  optically-pumped,  spin-polarized  atoms  in  the  presence 
of  optical  or  microwave/radio  frequency  radiation,  static  magnetic  fields,  and  various  types 
of  collision  processes.  The  book  shows  how  to  treat  real  atoms,  with  many  sublevels  in 
the  ground  state  and  optically  excited  state,  with  the  same  degree  of  rigor  as  the  idealized 
“two-level  atom,”  which  is  widely  discussed  because  of  its  simplicity,  but  which  is  of  limited 
use  for  modeling  quantitative  details  for  real  atoms. 

A  particularly  useful  way  to  describe  spin-polarized  atoms  is  with  a  density  matrix  p. 
The  book  explains  how  the  density  matrix  can  be  represented  by  a  column  vector  | p)  in 
Liouville  space,  with  a  time-evolution  equation 

>  =  -C|p).  (1) 

Often  the  the  evolution  matrix  G  is  independent  of  | p).  For  systems  with  oscillation 
perturbations  like  the  effects  of  microwave  radiation  or  laser  light,  a  generalized  “rotating 
coordinate  system”  can  often  be  found  for  which  G  can  be  taken  to  be  time  independent. 
So  (1),  often  called  the  “master  equation”  is  often  a  first-order,  vector  linear  equation  with 
time-independent  matrix  G  of  coefficients.  The  challenge  is  that  the  column  vector  | p) 
representing  the  density-matrix  p  has  a  very  large  dimension,  256  in  the  case  of  ground- 
state  133Cs  atoms,  and  the  dimension  can  easily  exceed  1000  if  elements  of  the  excited  states 
are  included.  As  stressed  in  the  book,  modern  mathematical  software  like  MATLAB  is 
well  adapted  to  calculations  in  such  high-dimensional  spaces. 

Unlike  the  analogous,  Schroedinger  equation 

jtW  =  -ihHW).  (2) 

which  evolves  under  the  influence  of  Hamiltonian  matrix  H  which  is  Hermitian,  //  I  =  II . 
the  damping  operator  G  of  the  master  equation  (1)  is  neither  Hermitian  nor  “normal,” 
that  is 

and  GG^G^G.  (3) 

We  assume  that  the  evolution  matrix  G  of  (1)  has  right  (column  or  ket)  eigenvectors  | g) 
and  eigenvalues  g  defined  by 

G\g)  =  g\g).  (4) 

Since  G  is  not  Hermitian  the  eigenvalues  g  of  G  need  not  be  real, 

9^9*-  (5) 

and  the  left  (row  or  bra)  eigenvectors,  which  we  denote  by  a  double  left  parenthesis  as  ((g\ 
and  define  by 

(g\G  =  {g\g>  with  {g\  ±  |#)t.  (6) 

need  not  be  Hermitian  conjugates  |^)t  of  the  right  the  eigenvector  | g),  defined  by  (4). 
The  situation  is  reminiscent  of  x-ray  crystallography,  where  it  is  customary  to  introduce 
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reciprocal  lattice  vectors  (left  eigenvectors,  row  eigenvectors)  that  differ  from  the  lattice 
vectors  (right  eigenvectors,  column  eigenvectors)  to  efficiently  analyze  diffraction  patterns. 

Systematically  building  on  the  simple  ideas  noted  above  the,  book  shows  how  to  use 
modern  mathematical  software  like  MATLAB  to  efficiently  solve  very  complicated  optical 
pumping  problems  where  the  density  matrix  has  hundreds  of  elements.  The  notation  of  the 
book  was  chosen  to  facilitate  computer  coding,  so  that  a  symbolic  equation  involving  vec¬ 
tors  or  matrices  looks  almost  the  same  when  converted  to  a  line  of  code  in  MATLAB.  The 
book  also  contains  many  software  packages  for  optical  pumping  that  can  be  downloaded 
from  a  permanent  website  at  the  University  of  Wisconsin  for  use  in  specialized  problems. 

Velocity  changing  collisions  of  spin-polarized  atoms.  One  of  the  significant 
results  of  AFOSR  research  was  the  principal  investigator’s  invention  of  the  the  “sodium 
guide  star,”  for  sensing  and  correcting  for  the  effects  of  optical  turbulence  on  propagation 
of  lasers,  for  sharpening  ground  based  imaging  of  space  objects,  and  for  sharpening  the 
resolution  of  large,  ground-based  astronomical  telescopes.  The  Na  guidestar  concept  was 
heavily  classified  throughout  the  1980’s,  when  extensive  experimental  work  led  by  Dr. 
Bob  Fugate  at  Kirtland  Air  Force  Base  demonstrated  its  usefulness.  At  altitudes  between 
about  90  and  100  km,  there  are  free  Na  atoms  in  sufficient  abundance  to  backscatter  a 
significant  fraction  of  laser  light  tuned  to  the  center  of  the  Na  resonance  line.  The  Na 
atoms  are  produced  by  the  vaporization  of  micro-meteoroids  that  constantly  bombard  the 
earth.  The  altitude  of  the  Na  layer  is  where  the  micro-meteoroids  get  hot  enough  from 
” reentry”  heating  to  volatilize  most  of  their  sodium  atoms  or  to  ’’burn  up”  entirely.  All 
large  modern  telescopes  now  include  a  sodium  guidestar  system,  which  consists  of  a  laser 
located  at  the  observatory  that  can  be  tuned  to  the  590  nm  D2  line  of  sodium  atoms  and 
directed  toward  the  star  or  galaxy  of  interest.  The  backscattered  laser  light  from  these 
atoms  can  serves  as  an  artificial  bright  star  to  probe  the  atmospheric  turbulence  and  close 
feedback  loops  to  correct  for  the  degradation  in  spatial  resolution  or  “  seeing”  that  is 
caused  by  the  turbulence. 

Although  Na  guidestar  systems  work  well  in  astronomical  applications,  there  are  no 
return  photons  to  spare.  One  can  use  physics  to  increase  the  backscatter  without  using 
more  powerful  lasers.  One  of  the  most  useful  ways  to  increase  Na  guidestar  signals  is 
to  optically  pump  atoms  into  spin  sublevels  with  large  backscatter  cross  sections.  The 
physics  related  to  this  pumping  differs  quite  substantially  from  that  encountered  in  most 
laboratory  experiments  on  optical  pumping,  which  are  normally  done  at  much  larger  gas 
pressures,  as  in  the  case  of  gas-cell  optical  pumping,  or  done  with  no  background  gas  at 
all,  as  in  the  case  of  optically  trapped,  cold  atoms.  At  the  altitude  of  the  sodium  layer, 
the  the  total  number  density  the  atmospheric  nitrogen  and  oxygen  molecules  and  the 
relatively  small  number  of  dissociated  oxygen  atoms  is  only  about  N  —  1014  cm-3.  The 
rate  coefficients  for  collisions  of  the  Na  atoms  on  the  residual  gas  molecules  and  atoms 
are  of  order  (va)  =  10“ 10  cm3  s_1,  so  collision  rate  is  only  7  =  N(va)  —  104  s_1,  much 
slower  than  the  natural  radiative  decay  rate  of  excited  Na  atoms,  1/r  =  6  x  10'  s-1  for 
the  radiative  lifetime  r  =  16  ns  of  the  first  excited  state  of  Na. 

Because  of  the  slow  collision  rate,  the  overall  width  of  the  optical  absorption  cross 
section  will  be  almost  unaffected  by  collision  broadening  and  will  be  almost  completely 
determined  by  hyperfine  structure,  natural  radiative  broadening  and  Doppler  broadening. 
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When  pumped  by  nearly  monochromatic  ground-based  lasers,  only  those  Na  atoms  moving 
with  “range  velocities”  v  that  Doppler  shift  them  into  resonance  will  absorb  laser  light.  To 
analyze  Na  guidestar  atoms,  one  must  therefore  keep  track  not  only  of  the  spin-state  of  the 
atoms,  but  also  of  the  velocity  state  [10] .  Denote  an  element  of  the  atomic  density  matrix 
for  atoms  in  the  ground  state  by  pr,  where  r  denotes  the  velocity-independent  “relaxation 
rate”  of  the  atom.  To  include  both  population  and  coherent  components  of  the  density 
matrix,  r  is  taken  to  be  a  complex  number  with  the  real  part  representing  relaxation  and 
the  imaginary  part  representing  the  spin-precession  frequency.  A  single  collision  can  cause 
a  large  change  in  velocity,  so  it  is  not  possible  to  describe  the  evolution  of  pr  with  a  partial 
differential  equation.  Instead  the  evolution  is  described  by  the  integral  equation  [10], 

q  r°° 

-pr(v,t)  =  -J  {y\Kr\v')pr{y' ,t)dv'  +  Sr{y).  (7) 


Here  {y\Kr\v')  is  the  kernel  of  the  integral  transform  that  describes  velocity-changing 
collisions,  and  Sr(v)  is  the  rate  of  creation  of  the  density  matrix  element  pr  with  relaxation 
rate  r  and  velocity  v.  For  simplicity,  we  have  taken  Sr  to  be  independent  of  time,  but 

it  is  straightforward  to  consider  time- dependent  rates.  The  dimensionless  velocity  v  is 
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measured  in  units  for  which  a  Maxwellian  distribution  would  be  e~v  . 

For  weak  collisions,  each  of  which  changes  the  velocity  by  a  small  amount,  the  integral 
equation  (7)  can  be  well  approximated  by  a  differential  equation, 


d_ 

dt 


Pr 


1  d2 
2diR 


d_ 

dv 


v  +  r 


Pr  "F  Sr> 


(8) 


The  differential  equation  (8)  represents  forced  diffusion  toward  a  steady-state,  Maxwellian 

2 

velocity  distribution,  e~v  ,  since 


Id2  d  \  _v2 

+  e  =0. 


2  dv 2  dv 

For  steady  state,  when  dpr/dt  =  0,  the  formal  solution  to  (8)  is 

/OO 

(v\Gr \v')Sr{v')dv' . 

-oo 

The  Green’s  function  for  the  velocity-diffusion  equation  (8)  is  the  solution  to 

1  d2  d 


(9) 


(10) 


2  dv2 


-^-v  +  r)  (v\Gr\v')  —  S(v  —  v'). 
dv  ) 


(11) 


In  an  important  advance,  Happer  and  Morgan  [10]  showed  that  one  can  write  the  Green’s 
function  defined  by  (11)  in  the  closed  form 

2  r  2 

(v\Gr\v')  =  ^^r(r)et;  Rr(—  v<)Rr(v>).  (12) 
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where  v<  is  the  lesser  of  v  and  v'  and  v>  is  the  greater.  T(r)  is  Euler’s  gamma  function  for 
the  complex  argument  r.  The  “right”  functions  Rr  are  defined,  except  for  normalization, 
by  the  differential  equation 


(13) 


with  the  boundary  condition 


Id  .  .  . 

- - v  Rr(v)  — >  0,  as  v  — >  oo. 

2  dv  1  w 


(14) 


The  right  function  is  entire  (that  is,  it  has  no  singularities  for  finite  values  of  v  in  the 
complex  v  plane).  A  convenient  expression  for  Rr.  is  the  power  series  [10] 


Rr(v)  =  ^ 


V^i-vy 


n!2r_nT( 


(15) 


n= 0 


Cusp  kernels.  There  have  been  very  few  laboratory  studies  of  optical  pumping  under 
conditions  like  those  of  optically  excited  Na  atoms  at  an  altitude  of  about  100  km.  The 
residual  atmospheric  pressure,  on  the  order  of  millitorr,  is  very  small  by  the  standards  of 
laboratory  gas  cell  experiments,  but  it  is  still  much  higher  than  the  pressures  for  work  with 
laser  cooled  and  trapped  atoms.  To  analyze  the  few  experiments  that  have  been  done,  one 
assumes  that  a  polarized  element  of  the  density  matrix,  p(v,t),  evolves  according  to  the 
integral  equation, 

d  t* 

—  p{v,t)  =  -('Tad  +  'YVd)p(v,t)  +7„d  /  W(v,v')p(v',t)dv'  +  P(v,t).  (16) 

The  first  term  on  the  right  of  (16)  represents  losses  of  the  atomic  polarization  at  at  velocity 
v  by  spin-damping  collisions  at  a  velocity  independent  rate  7 sd  and  losses  due  to  transfer  to 
other  velocity  groups  at  the  velocity-damping  rate  7 vd-  The  second  term  on  the  right  of  (16) 
represents  the  collisional  transfer  of  polarization  from  an  initial  velocity  v'  to  polarization 
observed  at  the  velocity  v.  The  third  term  represents  creation  of  polarization  by  optical 
pumping  at  the  rate  P(v,t).  Eq.  (16)  is  of  the  same  form  as  (7). 

The  most  widely  used  model  for  the  collision  kernel  W(v,v')  of  (16)  has  been  the 
Keilson-Storer  (KS)  kernel  [15] 


p-(v—av')2/b2 

wa(vy)  =  — —= — ,  (17) 

The  KS  kernel  has  a  single,  real  “memory  parameter”  a,  with  0  <  a  <  1,  and  a  corre¬ 
sponding  width  b  —  x/ 1  —  a2.  From  inspection  of  (17)  one  sees  that-  the  KS  kernel  assumes 
that  on  average,  a  single  collision  converts  atoms  with  a  velocity  v'  into  atoms  with  a 
Gaussian  distribution  of  velocities  v  and  a  variance  b2.  The  new  distribution  is  centered 
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on  the  velocity  av',  that  is,  the  atoms  “remember”  a  fraction  a  of  the  their  pre-collision 
velocity  v'. 

The  KS  distribution  appears  to  have  been  first  introduced  to  physics  in  1891  by  Raleigh 
[10]  in  connection  with  his  studies  of  velocity  distributions  produced  by  weak  collisions. 
It  turns  out  to  the  impulse  response  function  of  the  differential  equation  (8)  at  the  time 
r  =  —  in  a/r,  that  is,  it  is  the  distribution  of  velocities  v  at  time  r  >  0  that  would  have 
evolved  as  described  by  (8)  from  the  distribution  S(v  —  u)  at  initial  time  r  =  0. 

The  collision  kernel  W (v,  v')  of  (16)  is  the  probability  that  a  collison  transfers  an  atom 
with  an  initial  velocity  v'  to  a  final  velocity  between  v  and  v  +  dv.  It  must  therefore  satisfy 
the  normalization  condition 

POO 

/  dvW(v,v')  =  l.  (18) 


_  /2 

Collisions  should  also  transform  a  Maxwellian  distribution,  e  v  ,  of  initial  velocities  into 
the  same  Maxwellian  distribution  of  final  velocities,  so  a  physically  acceptable  kernel  must 
satisfy  the  constraint 

POO 


W(v,v')e~v'  dv'  =  e~v\ 


(19) 


The  KS  kernel  (17)  satisfies  (18)  and  (19),  and  it  has  the  advantage  of  a  simple  analytic 
form.  The  basic  problem  with  the  KS  kernel  is  that  it  does  not  resemble  kernels  deduced 
from  experimental  measurements.  When  plotted  against  the  final  velocity  v,  observed 
kernels  have  a  sharp  cusp  centered  on  the  initial  velocity  v’ .  They  do  not  look  at  all  like 
the  Gaussian  KS  kernel.  The  origin  of  the  cusp  is  the  large  number  of  weak  collisions  that 
correspond  to  large  impact  parameters. 

Working  with  students  Robert  Marsland,  Bart  McGuyer  and  Ben  Olsen,  Happer  de¬ 
veloped  a  powerful  new  “cusp  kernel”  [15]  that  looks  very  much  like  experimentally  de¬ 
termined  kernels,  satisfies  important  physical  constraints  (18)  and  (19),  and  has  a  simple 
analytic  form  determined  with  a  single  sharpness  parameter  s,  that  is  analogous  to  the 
memory  parameter  a  of  the  KS  kernel.  The  cusp  kernel  can  be  written  as  a  superposition 
of  KS  kernels  with  a  continuous  distribution  of  memory  parameters  a. 


Cs(v,v')  =  s  /  Wa(v,v')a8  1da. 


(20) 


For  large  sharpnesses  s>l,  the  cusp  kernel  is  dominated  by  realatively  sharp  KS  kernels 
of  memory  parameters  a  close  to  1.  It  turns  out  [15]  that  the  cusp  kernel  has  nearly  the 
same  functional  form  as  the  Green’s  function  (12), 


Cs(v,v') 


r(sy 


Rs(-n<)Rs(n>). 


(22) 


Many  other  advantages  of  cusp  kernels  are  described  in  reference  [15]  and  also  in  reference 
[17],  where  it  is  shown  that  cusp  kernels  give  excellent  fits  to  laboratory  measurements 
of  atoms  that  are  optically  pumped  under  conditions  similar  to  those  of  guidestar  sodium 
atoms. 
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